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Abstract 



For a system of two-chain spin ladders, the ground state for weak interlad- 
der coupling is the spin-liquid state of the isolated ladder, but is an ordered 
antiferromagnet (AF) for sufficiently large interactions. We generalize the 
bond-operator mean-field theory to describe both regimes, and to focus on 
the transition between them. In the AF phase near the quantum critical 
point (QCP) we find both spin waves and a low-lying but massive amplitude 
mode which is absent in a conventional AF. The static susceptibility has the 
form x(T) = Xo + aT 2 , with xo small for a system near criticality. We con- 
sider the dynamical properties to examine novel features due to the presence 
of the amplitude mode, and compute the dynamic structure factor. LaCu02.5 
is thought to be such an unconventional AF, whose ordered phase is located 
very close to the QCP of the transition to the spin liquid. From the Neel 
temperature we deduce the interladder coupling, the small ordered moment 
and the gap in the amplitude mode. The dynamical properties unique to 

near-critical AFs are expected to be observable in LaCu02.5- 
PACS numbers: 75.10.Jm, 75.30.Kz, 75.40.Cx, 75.50.Ee 
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I. INTRODUCTION 



The study of spin ladder systems |T| has made rapid progress since their synthesis 0] and 
identification. || Among the problems investigated both theoretically and experimentally 
are the nature of the ground state, which is dramatically different for ladders composed of 
odd and even numbers of chains, the behavior of holes in ladders, including their super- 
conductivity, ||] and the properties introduced by doping of spinless impurities. Further 
novel physical phenomena are found in coupled spin ladders, where the magnetic structure 
of even-chain ladders with unfrustrated coupling, is of particular interest because the ground 
state is expected to be a spin liquid at sufficiently weak coupling, but to order magnetically 
when the interladder interactions are strong. This system realizes the conditions required 
to investigate the QCP, introduced by Chakravarty et al. || in the context of a nonlin- 
ear cT-model description of the two-dimensional Heisenberg AF, and discussed extensively 
and more generally in Refs. and 0. In the language of these analyses, the interladder 
coupling is a parameter which tunes the system from the renormalized classical regime of 
long-ranged magnetic order, through the QCP to the one- dimensional limit where the spins 
are disordered by quantum fluctuations (quantum disordered). 

Such a structure is represented in three dimensions by the two-chain ladder compound 
LaCu02.5- H Initial experimental studies gave contradictory results on the nature of the 
ground state, as static susceptibility measurements |§ suggested a spin liquid state with 
a spin gap, while nuclear magnetic resonance (NMR) || and muon spin resonance (/iSR) 
i0| measurements indicated a transition to a magnetically ordered phase. In a brief study 



of both electronic and magnetic properties of the material, it was proposed |TI] that the 
conflicting observations could be reconciled if the system was located near the QCP of the 
transition between the two regimes, and this scenario was supported by detailed numerical 



simulations performed by Troyer et al. [12 



In this work we will analyze the properties of the system on both sides of the critical point, 
employing a mean-field approximation to the bond-operator technique which is generalized 
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to the magnetically ordered regime. From the evolution of the mode structure through 
the transition, we may compute both the thermodynamics and the dynamical magnetic 
properties of the system. In the disordered phase, the excitations are triply degenerate 
magnons with a spin gap which vanishes on approach to the QCP. In the AF ordered phase, 
we find that the modes evolve into two spin-wave excitations representing rotations of the 
ordered moment, accompanied by an amplitude mode corresponding to fluctuations in its 
magnitude. This latter mode has a gap which grows continuously with the moment, so it will 
be low-lying only near the QCP, and will develop far from the transition into a high-lying 
excitation which is ignored in the dynamics of a conventional AF. 

To characterize the ordered system close to the QCP we focus on the static susceptibility 
x(T), which vanishes at T = on approach to the transition. The results for x(T) in 
the mean-field theory agree well both qualitatively and quantitatively by highly accurate 
Quantum Monte Carlo studies carried out by Troyer et al. [[HJ Both sets of calculations 
support the proposal that LaCu0 2 5 is indeed an AF close to the QCP. 

We examine the dynamic structure factor to find features associated with the presence of 
an amplitude mode. This mode couples to neutrons in the same way as the transverse excita- 
tions, and we show that it should be observable in inelastic neutron scattering measurements 
on single crystals. Raman light scattering by the magnetic excitations is also discussed, but 
in this case it is difficult to find signals which may be ascribed unambiguously to this mode. 

The outline of this paper is as follows. In section II we develop the bond-operator 
formalism for application in the ordered magnetic regime, to describe the mode structure 
of the magnetic excitations and the zero-temperature phase transition. In section III we 
discuss briefly the statistics of the magnon excitations, in order to solve the mean-field 
equations at finite temperatures and to deduce the static susceptibility on both sides of the 
transition. We consider in section IV the dynamical magnetic properties of the system in 
the vicinity of the QCP, and calculate the dynamic structure factor for comparison with 
proposed experiments. Section V contains our conclusions and a brief discussion. 
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II. ORDERED MAGNETIC PHASE 



We begin by applying the bond-operator technique UT3|,[H| in the regime where the in 



terladder superexchange coupling is sufficiently strong to stabilize a magnetically ordered 
state. We will restrict ourselves to the case of simple AF order both along and perpendicular 
to the ladders. However, we shall return in section IV to a recent alternative proposition, 
and show that this gives rise only to notational differences. 
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Following Ref. p~3|| , the spin degrees of freedom on each ladder rung are represented by 
the bond operators 

^(Itl)-llt)), 
-^=(ITT)-iU}), 

^(in> + iu», (i) 

^(IU) + UT», 

where the arrows denote the direction of the left and right spins. The presence of a quadratic 
term with negative coefficient in the transformed Heisenberg Hamiltonian ensures that the 
singlet operators condense, so that one may take (sj) = s on each rung. AF ordering along 
the ladders occurs when, in addition, one of the triplet operators has a non-zero expectation 
value with alternating sign. This we take to be the component t z , which may be written as 

t iz = {-i)H+i iz , (2) 

where t z contains fluctuations about the average value t. The dynamics of these fluctuations 
may not be neglected, as in the treatment of the s degree of freedom, because their dispersion 
is strong. With this choice of basis states, the coexistence of a finite t with the condensed 
singlet s on a rung corresponds to an AF alternation along the ladder of states | ||) and 
| |t) °f fixed, staggered moment t, admixed with a singlet component of weight s — t. 

In the disordered or spin-liquid regime the three triplet modes remain degenerate and 
massive, a situation represented schematically in Fig. 1(a). The spin gap vanishes as the 



coupling approaches the QCP, where the three modes of the lower branch are spin-wave- 
like, with a linear dispersion about the wave vector of the magnetic ordering (Fig. 1(b)). 
In AF the ordered phase there is a spontaneous breaking of rotational symmetry, which by 
Goldstone's theorem is accompanied by massless excitations. These will be two spin waves, 
which represent rotations, or transverse oscillations, of the finite staggered moment. The 
third mode corresponds to fluctuations in the amplitude of the moment (longitudinal) and 
will acquire a finite gap, as shown in Fig. 1(c). In a conventional AF system one observes 
only the two spin waves, but close to the QCP the amplitude mode will also be of low energy 
at the zone center. 

The real-space axes, not to be confused with spin-space labels introduced above, are 
chosen with z along the ladder direction and x and y the directions of interladder cou- 



pling. Using the reduced unit cell [|12| in the x- and y- dimensions, but with doubling of the 
real-space structure in the z-direction to accommodate AF ordering along the ladder, the 
Hamiltonian of the three-dimensionally coupled system, 

H = J Sij.S r j + XJ S m j.S mj '_|_2 

j j,m=l,r 

(3) 

j 

after transformation to bond operators is written as 

H = H + H 1 + H 2 + H 3 + H A , (4) 
where the terms have the following origins. 



j,a m=l,2 



-EE Hm{sfsj + tjlq a -i) (5) 

j,a 771=1,2 

contains the contribution from the ladder rung interactions, and also the constraint which 
restricts the spin states on each rung to a singlet or one of three triplets. Here j is a unit 
cell index and m an index for the two types of rung in each cell. 
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Hi = i a j e e (cfe^SV + ciOr ^ m + + / + ^- c -) ( 6 ) 

j.a m=l,2 

and 

ri _ h j V \ " (i m \ j.m+lfj.m+1 t m j.mt + TO+lf j.ro+1 j.m , tt _ \ /y\ 
"2 — 4^ 2^ \ i j,a L j+z,f3 L j+z,a L j,P i j,a i j+z,a i j+z,f3 i j,(3~^ rl - c -) \< ) 

j,<x£P m=l,2 

are the terms quadratic and quartic in £ operators corresponding to ladder leg interactions; 
H 2 may not be discarded in the ordered system because it will contribute terms of order £ 
to the dynamics of the modes a = x,y. 

M — _IV \ " \ " (iM-f-m mf m 

j',q r]=±x,y m=l,2 

+ tj,atj+r],a S j+ri S j + H.C.J (8) 



and 



17 _ 1\' TV \ \ " ( + m\,m\ ,m ,m 
n 4 — s AJ2_^ 2^ 2^ \ l j,a l j+ V ,f3 l j+ V ,a l j,f3 
j,a rj=±x,y m=l,2 

_ f mUm\ f m f m , tj \ /n\ 
L j,a L j+ri,a l/ j+r),l3 b j,f3 ' n.U.) {V ) 

are the quadratic and quartic contributions from the interladder coupling terms, which link 
rungs of the same type m. 

In the mean-field approximation the singlet expectation value (sf 1 ) = s and Lagrange 
multiplier /x™ = fi are taken to be the same on all rungs (i, m), while the triplet expectation 
value (£™ ) = — (— l) m £ alternates. The constraint term incorporates a reduction of s due 
both to the presence of the £ term and to quadratic fluctuations of all three magnon modes. 
A direct coupling of the longitudinal (t z ) fluctuations to the magnitude of s is also relevant, 
and appears in the t* z t z s*s terms as 

l j,z b 3+z,Z°j°3 — b j : z L J+z,z y± c j 1 z L 3, z J \ L t/ j+z,z b 3+z,Z J 

- (t + P jiZ )(-t + t j+ z, Z ) (s~ ^Ji,z) (s-^+iJj+S*) (10) 
= — S 2 t + S 2 ij iZ tj + z,z + |* (fj,ztj,z + tj+z,ztj+z,z) + ^(t 4 ). 

The first line follows from substitution of the number operators into the constraint and the 
second from Eq. In the last line, the first term is the classical part, the second is a 
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fluctuation term which will appear in off-diagonal components of the Hamiltonian matrix 
for all leg and interrung coupling combinations in H and for all three polarizations t a , and 
the third is diagonal in the magnon operators so appears as a mass term unique to the t z 
modes. Note that the mass grows with the ordered moment t, as would be expected. 
The mean-field Hamiltonian in reciprocal space may be cast in the form 



H m (fi,s, t) = N {-\Js 2 + \Jt 2 - lis 2 - tit 2 + 11- 2JsV(A + A') 

+ E ( E [A£C f c + A£ (C f *-L + 0- k .)] (ii) 

ka [m=l,2 

+ [K^JL + K (#A a + tlj 2 - ka )] + [i - 2]} , 

in which N denotes the total number of ladder rungs, the coefficients for the transverse 
modes a — a = (x, y) are 

A£ = \j-n- ±\' J{s 2 -f){cosk x + cos k y ), (12) 
A£ = \ A' J(s 2 + t) (cos k x + cos ky), (13) 



h 1 * = J(s 2 -t)cos\k z (14) 



and 



A k = \ J{s 2 + t 2 )cos\k z , (15) 
and for the amplitude modes a — z, 

Al = \j-/i + 2Jt 2 (X + A') - \\' Js 2 (cos k x + cos k y ), (16) 

A k = \XJs 2 {cosk x + cosk y ), (17) 

and 

A' k 2 = 2A / k 2 = \Js 2 cos\k z . (18) 
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The part of H m ( pTTj) quadratic in the triplet operators is diagonalized by Bogoliubov trans- 
formation, and can be reexpressed in terms of the appropriate quasiparticle operators 7^ 
as 

H m (n, s,t) = N (-§ Js 2 + \Jt 2 - fis 2 -fit + fi- 2Js 2 f(X + A')) 

-JV(iJ-/i) - \N (\J - ii + 2Jt 2 {\ + A') (19) 

+EE<4(iiit + l)- (2°) 

In calculating the mode frequencies, all terms fourth order in s and t are found to cancel, 
giving the easily factorized results 



(i J - /J, - 2uJs 2 a u k ) (I J - M + 2z/Jl 2 <)] 1/2 (21) 



and 



k; 



"(J J - // + 2 Jt 2 (A + A')) (\J-fi + 2Jt 2 {\ + A') - 2v Js 2 a u k )] 1/2 , (22) 
where 

= A cos \k z ± |A'(cos fc x + cos fc y ). (23) 

The equation for a;^ 2 represents two branches for the amplitude mode in the doubled Brillouin 
zone; folding back of the zone in the k x and k y dimensions returns the four branches in the 
Brillouin zone of the real material. The equation for u>^ a contains four doubly-degenerate 
branches, of which the most interesting is the lowest-lying, uj^ a with k x and k y in the reduced 



Brillouin zone. This yields the two rotation modes, and from the factorized form (pl|) it is 
clear that the spin-wave condition is the same as the vanishing of the spin gap which gave 
the critical coupling for the transition from the disordered side, | J — /x = 2Js 2 (X + A'). 

That the massless modes described by the bond operators are indeed true spin waves 
can be shown in two ways. First, considering the limit of fully developed magnetic order, 
s = t = -7= and 



S 



L0 } 



kcr 



A + A' 



1/2 



2Js 2 J(\ + \>) \\k 2 z + \{kl + k. 



1/2 



(24) 



in the limit of small k. A textbook derivation []T5|| of the mode spectrum for the Hamiltonian 
in Eq. (|3D, for a spins of arbitrary magnitude S and with four spins per unit cell, returns 
exactly the contents of both lines in Eq. fl24|), with the condition S = s 2 = ~, as required. 
More generally, in the spin basis of Eq. (P it is straightforward to show that a spin wave, 
represented by a staggered, transverse component of S x on each rung, may be represented 
by the bond-operator states — \t x ) — \s), and similarly for S y , showing its equivalence to a 
t x (t y ) bond-operator mode in the presence of the singlet. 

The average singlet and triplet occupations may be represented by the reduced variables 

2Js 2 



d t 



2Jt 2 



It { It' 

i J-/i jJ-V 

in terms of which the mode frequencies (|21~D and (^) are given by 



(25) 



kcr 



(i J - n) ^(l-vd s al){l + vd t al) 



(26) 



and 



w£z = (\J - fi) ^(1 + d& + A')) (1 + dt(X + A') - ud a a v k ), 



(27) 



while the spin-wave condition, which is maintained everywhere in the ordered regime, be- 
comes 



d s = (X + A') . 



(28) 



By substitution into the self-consistency equation for the disordered solution, [|IT| the critical 
coupling A' c is given implicitly by 

-1/2 



1 

aTa 



A + A' 



n m (a; k ), 



(29) 
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where in the reduced unit cell there is only one degenerate magnon branch. In the remainder 
of this section we will consider the system at zero temperature, so that the magnon thermal 
occupation factors n m (uj£) are unity. The full mean-field equations are 

6\l 



i 

3 

2 2 

k 



iy l + d t (\ + \>)-±vd s al 

2 v 2^(1 + dt(X + A')) (1 + MX + A') - ud B ay k ) 



y 1 + \vd t a v v — \vd s al- 



ku 2yJ(l - vd s al){\ + vd t a v k ) 
) = = -|J-^-2Jt 2 (A + A') 



uJaW + dtjX + X')) (31) 



kv 



2yJ(l + dt(X + A')) (1 + dt(X + A') - vd s al) 



E 



isJa^(l + vdtdQ 



t 2^ (I - ud s al)(l + vd t al) 



and 



(^> = = \ J - /I - 2Js\X + A') - J(A + A') 



Of 



4 

l 



iy 2J{X + X')(l + d t (X + X')-ud s ai) 
2 it 2 J {I + dt{X + A')) (1 + d t {X + A') - vd s a v k ) 



+ E v ^ a \S)- ~ vd s a£) 



ku 2^j (1 — vd s a u k ){l + vd t a u k ) 
in which the first three terms on the right hand side of each expression constitute the classical 
equations, and the remaining terms the quantum corrections. The notation J2' denotes J2- 
Unconstrained solution of these equations for /i, s, and t at fixed A and A' > X' c gives a 
set of parameter values which are completely unrelated to those on the disordered side of 
the QCP, meaning a first-order transition. However, such parameters correspond to modes 
which are in general all massive, as the spin-wave condition is not obeyed, and have no 
physical meaning in the present problem of an ordered magnet. This discontinuity between 
ordered and disordered regimes in an approach combining classical mean-field terms with 
quantum corrections has been encountered previously in similar circumstances. |7|, [T6|JT7|| 
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By the Goldstone Theorem, any breaking of a continuous symmetry must be accompanied 
by the presence of zero-energy excitations, in this case spin waves which are massless in 
the long- wavelength limit. We proceed by enforcing the third mean-field equation ( |32| ) 
classically: this is precisely the spin- wave condition Eq. (^), and serves both to ensure that 
the rotation modes are massless at the zone center and to fix the value of one of the three 
original variables. 

The remaining two mean-field equations (^,^l|) correspond to the pair solved in the 
disordered phase, now generalized to include finite t. These may be combined to yield a 
single equation for dt as a function of 5\' = X' — A' c , which as d t — > returns the critical 
point equation (|29|) for A' c . Solution of this equation gives a weak first-order transition which 
may be traced in a linear expansion about A' c to the fact that the small dt term acts as a 
cutoff in the k summation which provides a logarithmic contribution in three dimensions. 
The linearized equation has the form 

d t (A + -B In d t ) = C5X', (33) 

which as shown in Fig. 2, describes reentrant behavior of the onset of magnetic order. 
Numerically, the value t (Fig. 2) where the solution at fixed A' becomes double-valued is 
0.048. The A' axis is however grossly expanded, and the reentrance occurs within the range 
of values of interladder coupling 0.118 < A' < 0.121. Because the first-order nature of the 
transition is extremely weak, we may use the formalism developed in this section to describe 
the features of the ordered magnetic system both qualitatively an on a quantitative basis. 
Illustration of the zero-temperature solution for the variation of the ordered moment with 
the interladder coupling is deferred to the following section, where it is presented together 
with the solutions at finite temperatures. 
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III. FINITE-TEMPERATURE SOLUTION AND STATIC SUSCEPTIBILITY 



Magnon Statistics 

To compute the properties of the coupled ladder system at all temperatures, it is nec- 
essary first to discuss the statistics of the magnon excitations, which are contained in the 
thermal occupation function n m (u>^). [O] The essential feature is that despite the bosonic 



commutation relations [13[] of the single magnon operators, these do not have conventional 
bosonic statistics because of the constraint (||) on their number, and this has a profound 
effect on the thermodynamics of the system at intermediate and high temperatures. An 
approximate approach which gives a good account of the effects of the constrained magnon 
ocupation in the case of the isolated ladder (disordered system) was introduced by Tsunet- 
sugu and coworkers in Ref. and the reader is referred to this work for a detailed dis- 



cussion. Here we summarize the results of applying this treatment in the ordered magnetic 
regime. 

The following discussion is simplified by the implicit assumption of the presence of an 
anisotropy field stabilizing one particular spatial direction for the staggered moment, a point 
to which we shall return in more detail below when considering the susceptibility. This 
results in a small gap in the spin-wave spectrum, preventing the easy axis from reorienting 
in an applied field and allowing for a reduction of the energy of one of the modes. In 
an iV-rung system (^N unit cells) with both rotation (a = a = S = ±1) and amplitude 
(a = z = S = 0) modes possible at the same wave vector k, and in the presence of a 
magnetic field h, the free energy per site is given by 

/ = ~ In {1 + **(/?) + 2 cosh(/?/i)z CT (/3)} , (34) 

where z a {(3) = Yl'k exp(— ^C0k ja ), is the partition function for each mode type. The magnon 
thermal occupation functions per mode in zero field are 

^ k «) = e -PuL + 3 ( 35 ) 
12 



in the disordered phase where u>-^ a = uj-^ z = uj^, and in the ordered phase 

= i + e _^ kz + 2e _^ kg • (36) 

That in the latter two equations the occupation function for a mode z (a) contains a de- 
nominator dependent on the energy of a mode a (z) encodes the effects of the constraint. 
At high temperatures, each of these functions approaches the limiting value of ~ per mode. 
The form of the occupation function required in the mean-field equations and denoted in 
Ref. as n m (c<Jka) is that including the zero-point term, i.e. the analog of the bosonic 



\ coth Q/3a; kQ! J, and is given by n m (u; ka ) = \ + n m (u; ka ) 



Mean-Field Solutions 

Considering first the approach to the transition from the disordered side, in Fig. 3 
is shown the evolution of the spin gap A, the minimum of the triply degenerate massive 
magnon dispersion, with interladder coupling at temperatures between and 0.5 J. The zero- 
temperature graph is that shown already in Ref. [jTI]] . Solution of the mean- field equations 
at finite temperatures is marginally more complicated because the initial system of two 
variables may not be reduced to a single one in the same manner. As the temperature is 
increased, the disordered phase becomes more robust as might be expected, and a larger 
interladder coupling is required to stabilize magnetic order. The spin gap in the isolated 
ladder is found also to increase with temperature, although this rise is somewhat more linear 
above T = 0.1J than would be given by the simplest possible formulation, A = ^Aq + T 2 . 

In the ordered phase, the characteristic parameter varying with the coupling constant 
is t, which gives the extent of staggered moment formation. This is shown in Fig. 4 for 
the same temperatures. The ordered moment rises abruptly at the transition, with no 
obvious indication at any temperature of the reentrant behavior discussed in the previous 
section, and then less steeply thereafter. However, the logarithmic evolution of the ordered 
moment given in Eq. ( |3"BD is manifest in a In- In plot, and precludes the extraction of a 
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power-law mean-field exponent from data within a range A A' ~ 0.15 of the transition. At 
high values of the interladder coupling, the ordered moment falls short of the limiting value 
t — -j= as A' — > 1, showing that in the bond-operator description of the isotropic limit there 
remains some admixture of rung singlets with the ordered spins, a consequence of quantum 
fluctuations. 

Turning now to the physical situation of a fixed interladder coupling constant at variable 
temperature, it is clear from Fig. 3 that any system will be disordered at high temperature. 
By solution of a pair of equations analogous to the zero-temperature critical-point equation 
(f29l), this phase boundary is found in the plane of T and A', giving the Neel temperature 
shown in Fig. 5. Because in this figure is measured in units of J ~ 1400X, it is evident 
that the real material LaCu02.5, with a Neel temperature || of 117K, is located extremely 
close to the QCP (in fact at A' = 0.127). We will use for illustration the interladder coupling 
value A' = 0.13, for which the dispersion relations have already been sketched in Fig. 1(c). 
The dispersion curves in the physical Brillouin zone are now shown quantitatively in Fig. 
6 for this parameter choice, at which the ordered moment has the value t = 0.14 and 
T N = 0.105J. 



The physical quantity most readily measurable which gives direct information about the 
magnetic state of a spin system is the static susceptibility x(T), and for this reason it has 
been the subject of extensive discussions fllj|l8l and measurements in a variety of ladder 
systems. [p|,|19| The susceptibility per site is 



Static Susceptibility 



X(T) 



d 2 f 



(37) 



dh 2 



h=0 



which from Eq. fl3"4] ) yields in the disordered phase [|Ti| 



X(T) = (3 



(38) 
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This is illustrated over a wide temperature range in Fig. 7(a), where it is clear that the 
interladder coupling has an effect only at low T, while the remainder of the graph is well 
described by the numerical results for the isolated ladder. |TB| The situation at low tem- 
peratures is displayed in the inset, which shows that as the critical coupling is approached 
x(T) rises more rapidly with temperature. The nature of this increase may be studied by 
examining the function 



9(T)=T 2 ^fj, (39) 



in which x' denotes the temperature derivative. For the a system with spin gap JT8[ A, 
x{T) = T~ a e~ l3A and g = A — a>T, while for a power-law dependence x(T) = aT a , as 
might be expected in a critical regime, g = aT. Fig. 7(b) shows this function at low-T for 
several values of the coupling constant A', demonstrating clearly the evolution of the system 
with interladder coupling from a spin liquid with the spin gap Ao — 0.5J and power a = \ 
prefactor of the isolated ladder, to a quantum critical phase in three dimensions, where 
a = 2. These results have a straightforward interpretation from the thermal excitations of 
a dilute gas of (noninteracting) triplet magnons whose dispersion is quadratic about a spin 
gap A, and becomes linear as A — > 0. At the QCP the susceptibility is due to excitations 
of spin waves, with temperature dependence determined only by system dimensionality. 

The calculation of the susceptibility in the ordered phase may proceed in one of two 
ways: for an ideal system with no interaction between the spins and the real-space or lattice 
coordinates, the ordered moment simply reorients in an applied field to be perpendicular, 
i.e. a spin-flop transition. For a system with an easy-axis anisotropy, the spin waves acquire 
a mass proportional to the effective anisotropy field h a , and this allows the moment direction 
to remain stable under application of a parallel field. For a review of both situations see 
Ref. [p(| . The susceptibility is computed in the zero-field limit, which may always be taken 
as an external field smaller than the intrinsic anisotropy energy, so we restrict our attention 
to the simplest case of the susceptibility contribution from normal mode excitations for the 
case of pinned moment direction in a parallel field. 
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In addition to the part due to excitation of normal modes, the susceptibility in the 
ordered phase contains a constant part due to the presence of the finite moment, and so 
may be represented as 

x(T) = Xo (T) + Xcxc (T). (40) 

Xcxc(T) is computed by analogy with the disordered phase from the free energy of the modes 
which may be excited, and is given by 

vll (T)=0 ^ (41) 

Because this part has its leading contributions from the spin-wave modes with dispersion 
cJko-, it can be expected always to have a quadratic T dependence at low temperatures. The 
component due to application of a perpendicular field is similar, but slightly more involved 
because the evolution of the mode o>kz with field is not linear for a finite gap, a point which 
will be illustrated below. The numerator contains also a term in z z (j3), which is much smaller 
than the spin-wave contribution when the mass of the z mode is finite, and higher-order terms 
which cancel at the transition. The dominant feature of the temperature dependence remains 
the T 2 part due to the spin-wave contribution, with a smaller prefactor. On averaging over 
the crystallite directions in a polycrystalline sample, Xexc(^) = |xLc(^) + f/dicCO- 

Similarly, the static part in the anisotropy-pinned case may be written as Xo(T) = |xo 
for a polycrystalline sample, as only the transverse part has a finite value and this will be 
observed as an average over all crystallite orientations. \o± ma y De computed in the bond- 
operator formalism by considering the additional term in the mean-field Hamiltonian in the 
presence of a finite magnetic field h, 

E h- (s£ + Sl % + Sl + Sl) = -ze Q/37 £ (iX + t%t%) . (42) 

i i 

The evolution of the bond-operator modes of a ladder system with application of a magnetic 
field will be presented elsewhere. Here we state the result that in general a field component 
h a acts to create an ordered moment t m in the operators tp and t 7 , where a ^ /? ^ 7. In 
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the presence of a finite, staggered (t z ) = t, application of a small field h x stabilizes a finite, 
staggered ordered moment (—it y ) = t m . One of the spin-wave modes is unaffected by the field 
and the induced moment and one increases linearly with h, while the amplitude mode energy 
rises only quadratically with h. We find that t m oc th, and that the magnetization per site 
stabilized by a small field h x is m x = —^tmt, giving for the transverse part Xo = ^ t ■ 
As a function of interladder coupling, this contribution appears very similar to the the results 
in Fig. 5 for t. Xo vanishes approximately in the manner of a second-order transition as the 
temperature is increased towards T^. When the same formalism is applied to the case of a 
longitudinal field, the ordered moment is O (t^j and it is clear that xl is vanishing, so that 
the total x h &s T 2 contributions only 

In Fig. 8 is shown the full static susceptibility x(T) ( f40D for a variety of values of 
interladder coupling in the ordered phase, A' > A' c . The curves begin at the finite value of 
Xo(T = 0) calculated in the preceding paragraphs, and their variation with temperature is 
approximately quadratic because Xo(T) is nearly constant at low T, particularly for values of 
the interladder coupling not in close proximity to the critical point. Because the scale of xo is 
significantly smaller than that of Xexc, the second-order transition in Xo{T) at Tn appears as a 
minor downward cusp. The results compare very well with highly accurate Quantum Monte 
Carlo studies on very large systems by Troyer et al. [12], which show the transition from 



spin liquid to ordered magnet occurring very close to A' = 0.12, and continuous evolution of 
the susceptibility from an exponential form in the disordered phase to a quadratic variation 
at the critical point with an additive constant part which grows continuously on moving 
to the ordered side. Quantitatively, the magnitude of the peak susceptibility is some 15% 
greater in the mean-field theory, and that of the constant parts approximately 20% smaller 
at common values of the coupling in the ordered phase. One may fit the measured static 
susceptibility || to the form x(T) = a+bT 2 at low temperatures, and quadratic temperature 
dependence is found [|TJ] to give a good account of the data. The intrinsic susceptibility xo 
of the coupled ladder system is very difficult to extract from the constant a, as this requires 
detailed knowledge of core atomic susceptibility terms, and we note only the qualitative 
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result that all terms involved are very small, a further indication of the proximity of the 
system to the QCP. That the mean-field theory is in such generally good agreement with the 
numerical results, which for an unfrustrated spin system and with demonstrably negligible 
finite-size corrections can be taken to be essentially exact, is presumed to be primarily a 
consequence of the three-dimensionality of the ordered magnetic system. 



IV. DYNAMIC MAGNETIC PROPERTIES 

In the preceding section we have deduced that the LaCu02.5 system is magnetically 
ordered but located close to the QCP, as a result of which the transition temperature and 
ordered moment are small. One consequence is that the magnetic modes corresponding to 
fluctations in the amplitude of the ordered moment, which in a conventional magnet are 
high-lying and play no role in determining the low-energy properties of the system, should 
have only a small mass here. In this section we study the dynamic magnetic properties 
of such a system to isolate those features due to the presence of the amplitude mode and 
predict the ways in which these may be identified definitively by experiment. 

We focus directly on the dynamic structure factor measured by neutron scattering, which 
may be written following Ref. |21j in the form 



S*(q,w) = ^£(<y«p-g«&) £ i^'^*(q)Mq) 

ct/3 i,(ij),r,T' 



X 



dte-^(exp (-iq.rv(O)) exp (iq.r i+I . i . )T ,(t)))(^ +T (0)S'f +ri . +T ,(t)} ) (43) 



Here i and j denote unit cells, and r and r' are vectors specifying the locations of the atoms 
in each cell. F T (q) is the magnetic form factor which describes the spatial extent of the 
spin density around the site r, and g T is the Lande factor. In this expression the space and 
spin variables are partially factorized, and each expectation value can be separated into a 
constant part, corresponding to elastic scattering, and a time-dependent part responsible 
for inelastic processes. 

We consider first the elastic magnetic component, in order to isolate those Bragg peaks 



with the highest intensity about which one may measure the dynamical structure factor most 



readily, and also to discuss the nature of the magnetically ordered state. In Ref. ]TT| it was 
assumed that the interladder coupling in LaCu02.5 was AF. This assumption was made on 
the basis that ferromagnetic coupling, which relies on Hund's Rule coupling between orbitals 
on an atom in the superexchange path, is generally weaker than AF coupling energy scales, 
and that because the system is one of low symmetry is was likely that the latter would 
have appreciable components on certain paths. In fact this assumption was not tested in 
detail, and the tight-binding fit to the Local Density Approximation (LDA) bandstructure 
was made using a single Wannier orbital on each Cu site, which was taken to be a mixture 
(due to the low symmetry) of Cu Zd x 2_ y 2- and 3d 32 2_ r 2-centred orbitals in the pyramidal 
CuC-5 system. The good agreement obtained with this procedure appeared to justify the 
assumptions made, but did not rule out alternative parameter combinations. The results of 
the preceding section of this work provide a more reliable means of estimating the magnitude 
of the interladder coupling. Recently a quantitative effort has been made by Mizokawa et 
al. [^] to understand the superexchange parameters in LaCu0 2 .5- These authors deduce 
some of the parameters of the electronic structure from measurements made by Cu 2p 
core-level spectroscopy, and use these in Hartree-Fock calculations. They conclude that the 
magnetically ordered state of lowest energy is that where the AF coupled spins in each ladder 
are ferromagnetically coupled between the ladders, and that the magnitude of the interladder 
interaction is less than 0.1 J. Superexchange estimates involving competing exchange paths 
are difficult, and can involve significant errors due to computing small differences between 
large numbers. Nevertheless this is the most systematic study performed on the system to 
date. We note that the results of the preceding sections are essentially independent of the 
sign of the coupling ratio A', which acts only to exchange the branch indices v = ±. 

The elastic magnetic scattering is described by the static structure factor, the time- 
independent part of Eq. fl43|), and has finite components 

S*(G) = N(Sr (l - (G.z)L) [f<7F(q)]V(G) (44) 
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only at the reciprocal lattice vectors G, which are the magnetic Bragg peaks. Here z is 
the direction of the ordered spin moment in real space, and is not yet known. In cuprates, 
the coupling of the spin system to the lattice which determines this direction is generally 
very weak, although by comparison with the parent phases of tetragonal high-temperature 
superconducting materials the spins may be expected to align along either the rungs or 
chains of the ladders. (S z ) is the magnitude of the ordered moment, and for the coupling 
value A' = 0.13 is of magnitude V2t ~ 0.2. Thus the square of this quantity yields only a 
4% effect, rendering the elastic scattering rather weak close to the QCP. 



is the structure function, obtained by summing over the sites in the unit cell, in which the 
convention used with the site labels r is such as to ensure that each j and j spins introduce 
factors of the opposite sign, corresponding to the spin density distribution around each 
site. The appropriate site labeling for one plane of Cu atoms is shown in Fig. 9(a) for the 
simple AF structure function, which we denote as Type I AF, and in Fig. 9(b) for the case 
with ferromagnetic coupling between the AF ladders, which we denote as Type II AF. It is 
straightforward to calculate ^(G), and the results are expressed by quoting the reciprocal 
lattice vector G of the magnetic Brillouin zone as (h, k, I), with o and e denoting odd and 
even integers respectively. For the Type I configuration 



while for Type II the results are identical with [x <-> y]; the displacement vector components 
a x and a y are shown in Fig. 9. These quantities, normalized to unity, are shown in Fig. 10 
for the two configurations in any reciprocal-space plane of odd / and for values of h and k 
between and 10. It is evident that there are Bragg peaks of nearly maximal amplitude, 
and so the best points in the reciprocal zone around which to investigate the dispersion of 
the dynamical modes would be for example (2, 0, o) and (0, 1, o) in Type I, or (3, 0, o), 




(45) 



64 cos 2 G x a x sin 2 G y a y (e, e, o) , (o, o, o) 
64 sin 2 G x a x cos 2 G y a y (e, o, o) , (o, e, o) 



(46) 







/ even 
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(0, 4, o) and (5, 1, 6) in Type II. Because the two configurations differ significantly in the 
locations of strong Bragg peaks, it should be possible by diffractometry to determine the 
sign of the interladder superexchange. 

Turning to the dynamical structure factor as it might be observed (in a single-crystal 
sample) around one Bragg peak, the factor of interest is the spin part fdte itJt (S^(0)S^(t)), 
from which the spatial dependence on the crystal structure has been removed. In order 
to compute this quantity in terms of the bond-operator eigenmodes, we take first the rung 
combinations 

SU t ) = T,^[SUt)±S na (t)], (47) 

i 

and then combine the variables on the two rungs in the reduced, AF ordered unit cell to 
give 



k 

S£ = s (C + t\ a ) ± 8 (t\ a + t 2 \ a ) . (48) 
The final expression for the structure factor 

/oo 
dte^(SZ(0)SlUt)) (49) 

is a sum over the three spin indices (for unpolarized neutron scattering) and all sign combi- 
nations of S , and its evaluation is aided by extensive cancellation of terms. For scattering 
studies at low temperatures, where very few excited states are occupied and these only 
singly, we may from their bosonic commutation relations |L3| compute the thermal expecta- 



tion value by summation over bosonic Matsubara frequencies. The two types of component 
which emerge from Eq. (HSf) are 



S R -(q, u) = J2^(s 2 + t 2 ) (cosh 28° - sinh2^) [n° + Q(u)] S (u% - \u\) (50) 

a 

S R+ ( q ,u)= J2 {(cosh2(^ +q -^) + l)nfK + l)5(^ +q -u;f-^) 

+ \ (cosh2(^ +q - ei) - i) K + ©M] K + e(w)] s (< +q + <4 - m)} , (si) 
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in which the hyperbolic trigonometric functions are the coefficients of the Bogoliubov trans- 
formation, JO]] and the are thermal occupation functions. Errors due to the constraint 
on magnon occupations are neglible at low temperatures, and may in fact be taken into 
account fully by a more general formulation of the scattering expression, which gives dif- 
ferent thermal occupation functions n k . S R ~(q, uj) ( pOD appears with coefficient (~s 2 + £ 2 ) 
for contributions from the a components of the spins, due to the presence of the ordered 
moment in t z , and with coefficient s 2 otherwise. It has the simple interpretation of scat- 
tering processes involving magnon creation or destruction, and will be seen as a broadened 
5-function line along (q, u; qQ ,) corresponding to the magnon dispersion relations. S R+ (c{,i 



(pl|) describes magnon-magnon scattering (first line) and pair creation or destruction pro- 
cesses (second line), each involving magnons of a different type, and so will have (a, a) and 
(a, z) components which could in principle be distinguished by polarized neutron scattering. 

The dynamic structure factor S R ~ (q, uj) for magnon creation (uj > 0) is shown in Fig. 
11, and as stated above is found to consist of a series of peaks at the dispersion relations 
of all magnon branches. A small broadening is inserted by hand in the calculation. In the 
disordered system (Fig. 11(a)) there are four triply degenerate branches in the physical 
unit cell, and the lowest of these has a gap at the bottom of its band (q = 0), as required 
for a spin liquid. In the ordered system (Fig. 11(b)), this gap has vanished, and branches 
are split into their a and z components. In the lowest-lying branch, the intensity of the 
spin-wave peak diverges as q — > and the magnetic Bragg peak is approached, while the 
corresponding z mode appears at the energy u = 0.29J at the chosen value of interladder 
coupling. This energy splitting is expected to be resolvable, and so the amplitude mode 
should be detected as a neighboring peak, or at least as a significant shoulder to the elastic 
peak, and to each of the dispersion branches. The component S r+ (q,uj) for two-magnon 
scattering and pair creation has an intensity between one and two orders of magnitude lower 
than that from magnon creation, and so will be barely detectable. Because inelastic neutron 
scattering couples directly to the one-magnon process, it is the best technique to observe a 
clear signature of the low-lying amplitude mode. 
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We have considered also the expected form of the intensity I{oS) in a Raman light scat- 
tering experiment on the magnetically ordered but nearly critical system. Unambiguous 
spectroscopic evidence for the amplitude mode, as in the magnon creation signal expected 
in the neutron scattering cross section, would be provided by one-magnon processes with a 
significant weight. Following Fleury and Loudon, |23| it is clear that in this system neither 
the magnetic-dipole nor electric-dipole first-order interactions detailed by these authors can 
provide appreciable signals. Concluding for completeness with second-order light scattering, 
we have computed the ((Sj.Sj)(Sj.Sj)) correlation function which gives the intensity in a 
spin-only model. |24j We find a "3J" peak due to scattering of zone-boundary magnons in 
the ladder (±/c 2 ) direction, and that this peak is broadened in the presence of an amplitude 
mode split off from the spin wave band in an AF system close to the QCP. However, mindful 
of the facts that measured spectra are not fully explained by two-magnon scattering consid- 
erations, and that these are found to be anomalously broad in other cuprate materials, we 
conclude that Raman scattering is unlikely to be a suitable probe of the dynamics of the 
amplitude mode. 



V. CONCLUSION 

We have presented a theoretical description of the three-dimensionally coupled spin lad- 
der system which is realized in the insulating compound LaCu02.5- This spin configuration 
is such that for weak interladder magnetic coupling the ground state is the spin-liquid phase 
characteristic of the isolated ladder, while for larger coupling values an unfrustrated, mag- 
netically ordered state is expected. We employ a representation based on singlet and triplet 
bond operators on each rung of the ladders to obtain a uniform description of both phases 
within the same framework. 

Within the mean-field approach to the bond-operator formalism, solution at zero tem- 
perature gives the interladder coupling at the QCP as J' = 0.121J. The transition is found 
to be very weakly first order, and the increase of the ordered moment to be logarithmic in 
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the coupling ratio A' — A' c in the vicinity of the transition, with the mean-field exponent 
recovered at higher A'. Finite-temperature solutions are obtained by incorporating the con- 
straint on the triplet excited states into an effective magnon statistics, and show the increase 
with temperature of the disordered regime and the critical coupling, allowing the deduction 
of the Neel temperature. 

The real material has an ordering temperature extremely small on the scale of the su- 
perexchange parameter J within the ladders, implying that it is located very close to the 
QCP on the ordered side. In this case the magnon mode corresponding to amplitude fluc- 
tuations of the ordered moment has only a small gap at the bottom of its band, and will 
contribute to the low-energy dynamic and thermodynamic properties of the system. The 
constant part of the static susceptibility is proportional to the square of the ordered moment 
and is small, reconciling measured susceptibility data with the fact that the system is or- 
dered. The part of the static susceptibility due to thermal excitation of modes of the system 
increases quadratically with temperature for spin waves in three dimensions. The results 
of the mean-field theory are supported by comparison with the detailed Quantum Monte 
Carlo studies of Troyer et al. []12j] . We have demonstrated that the low- lying amplitude 
mode, which is unique to this type of AF, contributes to the dynamical response function. 
Inelastic neutron scattering is an especially suitable probe for investigating this mode as it 
may access one-magnon excitation processes, and we have presented an explicit calculation 
of its effects for comparison with experiment. 

In closing, for a system so close to the QCP we may speculate on ways of controlling 
the tuning parameter represented by the interladder coupling in order to pass through the 
transition. As discussed in the preceding section, the superexchange processes contributing 
to the interladder magnetic interaction are not well understood, but the coupling is known 
to be a decreasing function of the bond angle away from 180°. It may be possible by 
application of hydrostatic pressure, or better uniaxial pressure along the x- or y-axis of a 
single crystal, to cause an alteration of this bond angle significant enough to measure, at 
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least as a raising or lowering of Tjv as deduced from i by NMR. As an alternative to physical 



pressure it is possible |25j also to apply chemical pressure by substituting other trivalent 
atoms for La. In LaCu02.5 the bond angle is found to increase on substitution of Y for La, 
while it is decreased in equal amount by substitution of Nd. We await with interest more 
detailed results on the evolution of the magnetic state, which may provide evidence of the 
system being moved through the QCP. However, for a random distribution of substituents 
there will be a random distribution of local distortions, from which one may deduce only 
an average value for the bond angle and interladder coupling, so that interpretation within 
the above framework will be less transparent. Finally, as we have shown in the preceding 
sections, the changes in many bulk quantities across the quantum critical regime remain 
essentially smooth, and it will be necessary to adopt a criterion based on the appearance of 
the staggered moment or the associated massive mode for unambiguous identification of the 
transition point. 
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FIGURES 

FIG. 1. Schematic representation of triplet excitation mode structure L0^ a in LaCuO"2.5 system, 
(a) Disordered (spin liquid) phase at low interladder coupling A' < A' c : triply-degenerate, massive 
modes, (b) Quantum Critical Point A' = A^,: triply-degenerate spin waves, (c) Magnetically 
Ordered phase at larger interladder couplings A' > A' c , with halved Brillouin zone. 

FIG. 2. Graphical representation of weak first-order part in transition of t(X'), from solution 
ofEq. (H). 

FIG. 3. Spin gap as a function of interladder coupling at several values of temperature, illus- 
trating monotonic decrease of A with increasing A' towards a vanishing at the point of transition 
from disordered to magnetically ordered phase. 

FIG. 4. Magnitude of ordered monent t as a function of interladder coupling at several values 
of temperature, illustrating abrupt increase from the weak first-order transition at A^., followed by 
more gradual progress towards nearly-full polarization. 

FIG. 5. Neel Temperature T/v as a function of interladder coupling. For LaCu02.5, the mea- 
sured value T/v = 117K < 0.1J || indicates very close proximity to the QCP. 

FIG. 6. (a) Static susceptibility x(T) for several values of interladder coupling in the disordered 
phase. Deviations from the behavior of the isolated ladder are observed only at temperatures 
T < X'J below the scale of the coupling energy. This regime is expanded in the inset, (b) The 
function f(T) ( |39| ) for four values of A' spanning the disordered phase, illustrating the change from 
exponentially activated to quadratic susceptibility with decreasing spin gap. 

FIG. 7. Magnon dispersion relations in physical Brillouin zone along the lines XT and TZ for 
magnetically ordered system A' = 0.13. All branches a are doubly degenerate, and the lower 
line is the amplitude mode with a small mass (0.29J here) near the QCP. 
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FIG. 8. Static susceptibility x(T) for several values of interladder coupling in the ordered phase. 
The constant part xo rises, at first rapidly, with A' away from A^.. Small downward cusps appear at 
T = Tjy, marking the second-order vanishing of xo- Isolated-ladder character is recovered at high 
temperature. 

FIG. 9. Unit cell structure and labeling for one plane of atoms in magnetically ordered material, 
for computation of structure function J-{G) (f45[). Intraladder coupling (perpendicular to plane of 
diagram) is always AF, whence unit-cell doubling, (a) Type I AF configuration, (b) Type II AF 
configuration. 

FIG. 10. Structure function J-(G) at reciprocal lattice points (h, k, Z=odd), for < h, k < 10. 
(a) Type I AF. (b) Type II AF. The full static structure factor requires also multiplication by the 
square of the ordered moment, and so will be rather small for a system close to the QCP. 

FIG. 11. Dynamic structure factor S R (q z ,uj) expected from inelastic neutron scattering, illus- 
trating contributions from magnon creation processes, which appear as a broadened <5-function at 
the dispersion branches uj a (q z ) (a) Disordered system, A' = 0.05: all four branches in physical 
Brillouin zone are triply degenerate, and there is a gap in excitation energy u to the first mode 
at q = 0. (b) Ordered system close to the critical point, A' = 0.13: branches are split into a and 
(massive) z modes. For the lowest branch, the spin- wave scattering amplitude diverges as q 2 as 
q — > 0, forming the elastic peak. The signature of the massive but low- lying amplitude mode is the 
peak dispersing from cu = 0.29J at q = 0. 
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